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We match the density of energy eigenstates of a local field theory with that of a random Hamil- 
tonian order by order in a Taylor expansion. Our earlier work showed the match to be sufficiently 
close to account for the success of field theory in describing our world, even if the underlying Hamil- 
tonian were simply a random one. In that work the assumed Lorentz symmetry of the field theory 
entered through the dispersion relation. Here we extend that work to consider a generalized disper- 
sion relation and show that the Lorentz symmetric case is preferred, in that the Lorentz symmetric 
dispersion relation give a better approximation to a random Hamiltonian than the other local dis- 
persion relations we considered. We concluded that a model of the physical world based on the 
approximate interpretation of a fundamental random Hamiltonian as a local field theory predicts 
Lorentz symmetry for the local physics. 



1. INTRODUCTION 

We consider the question "can a random Hamiltonian 
be interpreted as a field theory to a sufficiently good ap- 
proximation to describe the observed physical world?" . 
We posed this question in earlier work [l], Q from the 
point of view of a comparison of the energy eigenvalue 
spectrum, and found an affirmative answer. In that work 
we considered only field theories with a dispersion rela- 
tion given by Lorentz symmetry. The dispersion relation 
is the only way the presumed Lorentz invariance of the 
field theory entered into our calculations. In this work we 
extend our previous analysis to include a one parameter 
family of possible dispersion relations. We also include 
the dimension of space d (as was already analyzed in 0). 
Our analysis shows that choosing a Lorentz invariant dis- 
persion relation makes it easier to approximate the field 
theory density of states by that of a random Hamilto- 
nian. We thus conclude that in a physical picture where 
local field theory is just an approximate interpretation of 
a physical world which is fundamentally described by a 
random Hamiltonian, Lorentz invariant field theories are 
favored. 

Our past work on this topic was motivated by the 
"clock ambiguity" in quantum gravity. That motivation 
and how it takes us to the comparison of energy spectra 
is discussed at length in Q and 0. In this paper we 
focus more narrowly on extending our formalism to gen- 
eralized dispersion relations and interpreting our results. 
We feel that this technical work will be of interest to a 
wider group that may bring different motivations than 
our own. In section [5] we lay out our basic technical ap- 
proach in a manner that is fairly self-contained. The key 
new technical ingredient of this paper is our derivation of 
the density of states for field theories with a generalized 
dispersion relation. We present two such derivations, one 
using statistical mechanics (section[3]) and one using scal- 
ing arguments (section 2]) , and discuss the relationship 



between these two approaches in section [5l In section 
IH1 we integrate the generalized dispersion relation results 
into the framework presented in section [5] and show how 
this leads to our main conclusions. Section [7] reviews our 
conclusions and explores some interesting directions for 
further study. 



2. OUR BASIC APPROACH 

If the evolution of the world around us were most fun- 
damentally described by a random Hamiltonian, could 
that Hamiltonian be approximated to a sufficient degree 
by a field theory to allow us to interpret the physical 
world in terms of a field theory? In [l[ we addressed that 
question at the level of the spectrum of energy eigenval- 
ues. We found that the spectrum of a sufficiently large 
random Hamiltonian could approximate that of a free 
field theory over a sufficiently large range of energies to 
allow the former to be approximated by the latter on 
scales relevant to our data about the physical world. In 
this section we review that work with an eye toward the 
topic of this paper: The extension of that analysis to field 
theories with a generalized dispersion relation. 

The distribution of eigenvalues for a random Hamil- 
tonian, represented as an Nh x Nh Hermitian matrix, 
follows under quite general assumptions the Wigner semi- 
circle rule in the large Nh limit, for example, the distri- 
bution of eigenvalues of a large Hermitian matrix with 
elements drawn from a Gaussian distribution. 

On the other hand, the density of states for a free 
field theory grows at large energies like an exponential 
of a power of the energy. On the face of it, these two 
forms for dN/dE are dramatically different. In order 
to press forward with the comparison we introduced a 
general parametrization for the random Hamiltonian and 
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field theory spectral densities respectively: 
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The parameters Em and Es represent the maximum 
eigenvalue of the random Hamiltonian and an offset en- 
ergy between the two descriptions, fey (= 2n/L) is the 
resolution in fc-space set by putting the field theory in a 
box of size L and B, b, a and 7 are dimensionless param- 
eters. Since we have not yet explored the emergence of 
gravity in this picture we initially allow the energy off- 
set Eg to be a free parameter. The value of fey sets a 
scale at which continuum field theory breaks down. The 
standard expression 



N = exp S 



leads to 
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Which relates parameters in Eqn. [5] to the entropy S. 
Throughout this paper we we use units where K = c = 
k B = l- 

We expand both Eqns.[TJand[2]in a Taylor series around 
a given central energy £0 = puVu — 10 80 GeV~ (using 
the energy density pxj and volume Vu of the observed 
Universe today). We explore the discrepancies at each 
order in (E — Eg) to find the level of agreement between 
the two descriptions. 

Requiring exact equality at zeroth order sets the size of 
the space of the random Hamiltonian to be exponentially 
large: 
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More generally, this expression should be seen to only 
give a lower bound on Nh, since we currently only know 
upper bounds on kv (which gives a scale of breakdown 
of continuum field theory) . 

Requiring equality at first order (as well as equality 
at zeroth order) sets the offset energy Eg in terms of 
the energy of the Universe Eq by the following implicit 
expression: 
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Figure[T]gives examples of the two density of states curves 
equated at zeroth and first order according to the above 
analysis. Assuming equality and 0th and 1st order, the 
relative difference at second order is fixed and given by 
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FIG. 1: A plot of the density of energy eigenvalues for a 
random Hamiltonian using Eqn. [1] and a field theory using 
Eqn. [2] matching the zeroth and first order terms in a Taylor 
expansion around Eq (the vertical line). 



In [l| we considered specific values for the parameters 
in Eqn. [7] based on properties of the observed physical 
world and found that A2 is small for realistic values. We 
concluded that the hypothesis that the laws of physics 
could most fundamentally be given by a random Hamil- 
tonian, and just interpreted to a good approximation as 
a field theory had passed our "energy spectrum test" . 

The goal of this paper is to repeat the above analysis 
with a generalized dispersion relation used to get the field 
theory density of states. Since Lorentz invariance enters 
only though the dispersion relation, we hope to learn to 
what extent Lorentz invariant field theories are preferred 
when a random Hamiltonian is approximated by a field 
theory. To this end we introduce a general dispersion 
relation with two parameters, m and 9, of the following 
form: 



8{p)=p 



(8) 



Here £ represents the energy of a single particle with 
momentum p, and should not be confused with E, the 
energy at which we evaluate the density of states N(E). 
In general E represents the energy of multi-particle states 
(generally comprised of particles with with more than 
one energy). When the power g vanishes we recover the 
dispersion relation of relativistic massless particles and 
deviations from g — break Lorentz symmetry. The 
energy scale m is needed to keep the dimensions correct 
for nonzero values of g. A discussion of possible values 
of to, which is critical for the interpretation of this work, 
takes place in Section [6j 

There are a number of arguments to discard dispersion 
relations with negative values of g (see [3[ and references 
therein for related studies); a very direct one is as follows: 
If we allow g < 0, the group velocity, d£ /dp (proportional 
to p 9 ), becomes divergently large for packets with com- 
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ponents of low momentum. Thus, if we require signals 
to propagate at a finite speed (below a certain cutoff of 
accessible momenta) we are restricted to consider non- 
negative values of g. 

We conclude this section with some comments about 
additional assumptions and points of view implied by 
our framework. The formulas for the field theory den- 
sities of states quoted here (and generalized in the rest of 
this paper using Eqn. [8]) assume flat space (Minkowski 
spacetime for g=0). Following the approach developed 
in earlier work [l], Q and similar in spirit to what we 
have discussed here for field theory, we consider the pos- 
sibility that general relativity gives only an approximate 
description of physics that need only be valid around the 
observed state of the Universe. Thus the field theory 
in Minkowski space and its extensions which we consider 
here, without any explicit account made for gravity, seem 
like a reasonable place to start. To the extent that we can 
argue here that Lorentz invariant theories are favored we 
will be led to consider fields in different representations 
of the Lorentz group, including spin- 2. This can lead to 
the emergence of general relativity, as discussed briefly 
in section [71 

3. THE DENSITY OF STATES FOR A 
GENERALIZED DISPERSION RELATION: STAT 
MECH APPROACH 

3.1. Outline 



Throughout this paper we take degeneracy n = 1 for 
simplicity, without impacting our main points. For a 
free field theory in states of no net charge (and thus zero 
chemical potential, characteristic of the observed Uni- 
verse) the momentum distribution f(p) is given by 

= eW±T (13) 

where the upper/lower signs are for Fermi/Bose statistics 
respectively. 

Using Eqn. [5] in Eqn. Q2] gives 

= (TV* I 0FF ddp 

{2ir) d J e p(£J / T ± 1 
= T{Tm g ) d,{9+1) F(g,d) . (14) 
Here F is O(l) and is given by 

F{g,d) = -A-; / — - T d d z 

vy ; (2ir) d J e z9+1 ± 1 

To arrive at Eqn. [T51we changed variables in Eqn. Q3] us- 
ing z = p(Tm 9 ) 1 ^ 9+1 ' and then y = z 9+1 . We show 
explicitly in section [6] that F(g,d) is sufficiently slowly 
varying and close to unity that it can be safely neglected 
when building intuition about the expressions we are cal- 
culating. (Still, for completeness we keep it in for our 
final analysis.) 



For a thermodynamic system at temperature T in a 
volume V, the density of states is given by 



dN{E) 
dE 



dE 



where 



S = Vs = V 



p + P 
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(9) 
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gives the (extensive) entropy S in terms of the intensive 
quantities pressure P and energy density p. In this sec- 
tion we construct statistical mechanical expressions for 
S(T) and E(T) and then invert the latter to allow us to 
write 



S(E) = S(T(E)) 



(11) 



Although such a manipulation is only valid in the thermo- 
dynamic limit, we expect this to be a good approximation 
for a system the size of the observable Universe. 



3.2. Energy density 

For a homogeneous statistical system of degrees of free- 
dom of degeneracy n in d space dimensions the energy 
density is given by 



P(T) 



(27T) 



£(p)f(p)d d p . 



(12) 



3.3. Pressure 

There are some additional subtleties in defining pres- 
sure for a generalized dispersion relation. We start by 
rewriting the thermodynamic relation P = —dE/dV in 
the following way 

dE n f dp d£(p) 

W = ~ (2ir) d J W 



P 



W j 3v dp il6) 

Next, as in standard kinetic theory, we associate a vol- 
ume expansion SV — ASx to the pressure on the area 
A due to elastic reflection of a particle with momentum 
p x , that results in 8p x — —2p x A5x. For a expansion in d 
dimensions this generalizes to 

dp 2 



dV 

since Sx ~ Sr/y/d and p a 
P 2 n 



~d P > 
p/Vd. Thus, 

P^f(p)d d P . 
dp 



d(2n) d 

Using the dispersion relation of Eqn. [SJ namely, 
d£(p) 



dp 



we find 



(s + i)£(p) , 



2g + 2 



(17) 
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3.4. The statistical result 



4.2. The scaling result 



We find it convenient to write the volume V in terms 
of the inverse length given by 



kv = V 



-1/d 



(21) 



We are now in a position to complete the approach out- 
lined in section 13.11 using the expressions for P(T) and 
p(T) derived in the previous two sections. The result is 



S(E,g,d) 
d+2g 



2 9 + 1 
F d + g + l 



E 
kv 
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.(22) 



Note that both of the exponentiated factors in Eqn. [22] 
are greater than unity. The quantity kv sets the energy 
scale for the breakdown of field theory, and it should 
certainly be smaller than to for the generalized dispersion 
relation to have any meaning, leading to m/fcy > 1- The 
factor with E is discussed explicitly in section [2] 



Rescaling volume and entropy by A wc get from Eqn. 
and Eqn. [25] 
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(26) 



that is consistent only if the quantity in square brackets 
is independent of S. This implies 



S{E,g,d) = S (g,d) 



F \ d + g + l 



(27) 



The similarity with Eqn. [22] is suggestive of a close rela- 
tion to the statistical approach that will be explained in 
the next sectionlla. 



4. THE DENSITY OF STATES FOR A 
GENERALIZED DISPERSION RELATION: 
SCALING APPROACH 

4.1. Outline 

In this section we arrive at a result equivalent to 
Eqn. [22] via an alternative approach. We consider a ther- 
modynamic system characterized by its total energy, en- 
tropy, volume and energy scale to. We then make two 
assumptions: 

1. The energy is an extensive function of volume and 
entropy 

i.e., under a scaling of V and S by a factor A, 

E(XV, XS) = X E(V, S) . (23) 

2. The total energy is parametrized by the following 
ansatz: 

(24, 

where E ns stands for the energy in the case where 
there are no other scales (such as to) entering in the 
description of the system. This no-scale energy, 
by dimensional analysis, must be proportional to 
V- 1 ' 4 , i.e., 

V^ d E ns (V, S) = h(S) , (25) 

is a function of entropy only. 

Note that in Eqn. [M] we parametrized the dependence on 
the extra scale to as a power law in analogy with Eqn. [8] 



5. COMPARISON OF THE TWO APPROACHES 

The two approaches to obtain the density of states are 
designed to reflect different physical properties. The first 
one deals with variations of the energy density and pres- 
sure due to an altered dispersion relation in a way proper 
to statistical mechanics; using the representation of those 
quantities in terms of momentum distribution functions. 
The change in the dispersion relation is taken as the mea- 
sure of breaking of Lorentz invariance. The second ap- 
proach, in contrast, is centered in the response of the 
system to changes of scale, that is, in the breaking of 
conformality in the system. The results show, however a 
tight interplay of the two approaches. Indeed, if we con- 
sider a single particle state, we obtain a particular case 
of Eqn. [24] namely, 

E = £ ns {^Y , (28) 

where E ns is the no-scale energy for a single particle. But, 
since we assumed the total energy to be extensive, we can 
sum all single particle states (with their corresponding 
density, for example, in momentum space) to obtain the 
total energy, namely, 

E = Vj^- d J f(p)£{ P ) d d p , (29) 

where £(p) should be replaced with the expression on the 
rhs of Eqn. [2H1 Thus, by making the identifications 

Snsip) = P , (30) 

9 = 9, (31) 

we see that Eqn. [5S] coincides with the dispersion rela- 
tion Eqn. [8] and the energy density becomes equivalent 
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to the one derived in the statistical approach. Further- 
more, setting the undetermined factor So(g, d) in Eqn.l2"Tl 

g+1 

to (d + 2g + 2)F d +3+i /d, makes the entropy identical to 
the statistical counterpart in Eqn. [22] 

We note though that this last step is done simply by 
construction, and is not a derivation. The price to pay 
for the simplicity and generality of the scaling argument 
is that the coefficient So remains undetermined, while 
its numerical value is explicit in the statistical mechanics 
approach. 



6. ANALYSIS 

We are now ready to take up our original goal: The 
comparison of the field theory density of states with the 
Wigner semicircle. We have generalized the field theory 
density of states to include the parameter g which deter- 
mines the dispersion relation according to Eqn. [8] Since 
our results preserve the general form taken in Eqn. [2] for 
the field theory density of states, we can continue to use 
the formalism reviewed in section [5J 

Specifically, comparing Eqn. (32] with Eqns. [5] and [3J 
leads to 
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where 
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(35) 



Recall that A2 measures the fractional discrepancy be- 
tween the field theory and Wigner densities of states 
when evaluated over an energy range AE centered at 
energy Eo. The main new ingredient here is that A2 
now has a dependence on the dispersion relation via the 
parameter g. We start with an analysis of this new g 
dependence and then place this in the context of a full 
analysis of all the parameters in A2, the rest of which 
have already been analyzed in [l| and Q. 

The smaller the value of A2 , the better the field theory 
can be approximated by the Wigner semicircle (and thus 
by a random Hamiltonian) . We therefor expect values of 
g that lead to smaller A2 to be strongly favored in our 
scheme, because these values will lead much more easily 
to a field theoretic interpretation of a random Hamilto- 
nian. 

The dominant contribution to the g dependence of A2 
comes from the S(Eq) 2 factor in Eqn. [Ml In turn, the 
dominant contributions to the g dependence of S(Eq) 
come from the appearance of g in the exponents of Eq and 
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FIG. 2: The exponents of the E factor (solid) m factor 
(dashed) in Eqn. I35I evaluated at d £ {1,2,3,4} (height of 
the curve increases with with increasing d). Note that the E 
exponent decreases with increasing g while the m exponent 
grows. 



to in Eqn. [35] Figure [5] illustrates that the exponent of 
the E factor in Eqn. [35] decreases with increasing g, while 
the exponent of the to factor grows with g. The values of 
g that lead to smaller values of A2 thus depend on which 
factor dominates. Since (as discussed near Eqn. |2"2")) both 
exponentiated factors in Eqn. [35] are always greater than 
unity, the question of which behavior dominates is related 
in a simple way to the relative sizes of to, Eo and fcy. A 
helpful quantity is 



Evaluating Eqn. 



(£o4) 1/(d+1) 



]at m = m c gives 



S (to = to c , g, d) 



d + 2g 



2 L_ / TO, 

-Pd+g+l [ 



(36) 



(37) 



Equation [37] shows that when to is set to the critical 
value to c the two exponentiated factors in Eqn. [35]can be 
collected into a single factor which has no g dependence, 
as illustrated in Fig. [3] 

Figure [3] also shows that values of the dispersion rela- 
tion scale to less than m c make S decrease as a function 
of g and thus favor divergently large values of g, whereas 
to > m c causes g — to be favored. Since large g cor- 
responds to a non-local field theory, we conclude that to 
make a local field theoretic approximation to a random 
Hamiltonian one must choose to > m c in the dispersion 
relation. Having made that choice, the Lorentz symmet- 
ric case, g — is favored, in that g = minimizes A2 and 
thus can be best approximated by the random Hamilto- 
nian. 

For completeness we now examine the dependence of 
A2 on g from contributions other than the dominant E 
and to factors in S. We define the "prefactor" P from 
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FIG. 3: The entropy S(E ,g) for d = 3 for m/m c = 0.01, 
0.1, 1, 10 and 100 in order of lowest to highest curves. We 
see that taking m = m c removes the ^-dependence of S. For 
m > m c , S increases with g, whereas S decreases for m < m c 
. (Here E = W 80 GeV and kv = W 25 GeV) 




FIG. 4: The value of the prefactor P (defined in Eqn. [3Sj) 
as a function of g. The three curves show d = 2, d = 3 and 
d = 4, lowest to highest. The fact that P(g,d) = 0(1) for a 
wide range of values is why we can neglect P when building 
intuition about the dependence of A2 on g. 

Eqns. [3! and [35] as 

P(g,d)= d + 2! ! + 2 a 2 F^ . (38) 
a 

which contains additional ^-dependence not accounted 
for by the exponents of the E and m factors in Eqn. 1351 
Figure H| shows explicitly that P is 0(1) for a range of 
parameters as claimed. Still, for m very close to m c the 
shape of P(g) can have an impact, as illustrated in Fig. [5] 

Since we have argued that the g dependence of A2 
drives us to the the Lorentz symmetric case already con- 
sidered in earlier work, our analysis of the dependence 
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FIG. 5: This figure is similar to Fig. [3] but shows values of m 
much closer to m c (m/m c — 0.9, 0.95, 1, 1.05 and 1.25). For 
m so close to m c the non-trivial shape of the prefactor shown 
in Fig. U becomes visible and can lead to local minimum 
somewhat away from g = 0. That shifts the value of m for 
which = is preferred to a value slightly higher than m c . 



of A2 on the other parameters carries over unchanged. 
In particular, in Q we considered the same d depen- 
dence used here and came to the intriguing conclusion 
that the d — 3 case put A 2 right at the edge of values 
allowed by current data. To achieve that one had to take 
fey Ri m 7 = 10 25 GeV", where m 7 is the upper bound 
on the photon mass. Since kv denotes the scale where 
continuum field theory breaks down (not necessarily the 
scale of any actual volume V), this seems to be a very 
reasonable choice. Note also that in Q we simply took 
b = 1 in Eqn. [5] Here we have derived an explicit for- 
mula for b (given by Eqn. [32]) , which indeed is O(l) for 
the (? = case. 

As already discussed, A2 quantifies the degree to which 
the density of states curves for a field theory and a ran- 
dom Hamiltonian can over lap over an energy range AE. 
Figure [S] gives an explicit illustration of what is going 
on, where one can see how the two curves match more 
closely for the g = case. In this figure, some differences 
between the curves from one panel to the next reflect the 
different impacts of enforcing equality at zeroth and first 
order (Eqns. [5] and [6]) for different values of g. 



7. DISCUSSION AND CONCLUSIONS 

We have calculated the density of states for a free field 
theory with an arbitrary power law dispersion relation. 
This allowed us to extend our earlier work examining the 
degree to which a random Hamiltonian could be approx- 
imated by a local field theory. Our conclusion is that 
choosing the Lorentz invariant £{j>) = p dispersion rela- 
tion reduces the discrepancy between the energy eigen- 
value spectrum of a local field theory and that of the ran- 
dom Hamiltonian, as compared with the other dispersion 
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FIG. 6: NJj 1 dN /dE for random Hamiltonians (dashed lines) 
and field theory (solid lines) matching the zeroth and first or- 
der terms in an expansion around Eo as a function of E/Eo 
for g — (upper panel), g — 1 (center) and g — 2 (lower). The 
Lorentz symmetric (g = 0) case gives the broadest overlap, 
illustrating the main conclusion of this paper. The plots cor- 
respond to d = 3, p = 2, 7 = 1/2, k = 10-3S M , E = 8E M 
and m — Em- The parameters Nh and Es were solved for as 
explained in section [2] 

relations considered. In this sense, Lorentz symmetry is 
favored by our analysis. This is the concrete result from 
our work. 

As discussed in 0, , quite a few steps still need to be 
taken to check the viability of the "random Hamiltonian" 
picture. The work presented here offers additional hints 
about specific areas for further investigation. Imposing 



Lorentz invariance on the held theory will be the most 
general way to ensure the preferred dispersion relation. 
One then would expect spin-1 and spin-2 representation 
to turn up in the interpretation of a random Hamiltonian 
as a held theory. According to general arguments 0- 
[Toj these fields would lead to gauge forces and general 
relativity respectively. If we assume that the random 
Hamiltonian is the correct underlying theory, this picture 
should provide a systematic and hopefully testable way 
in which local held theory, gauge theory and gravity will 
break down. 

Looking more specifically at the technical results in 
this paper, we chose m > m c (with m c defined in Eqn. 
I5rj)) to avoid favoring arbitrarily large powers of momen- 
tum p in the dispersion relation 

£(p)=p(^y (39) 

which would give a nonlocal held theory. This led to a 
preference for g = which seems to eliminate m from the 
picture altogether. However, m c may show up in some 
higher order way. The two values considered in [l[ and Q 
for ky are ky = (Vu) 1 ' 11 (where Vu is the Hubble volume) 
and ky — m 7 . Continuing to use E = puV\j — lO so GeV 
as we have here and in our earlier work, these give 
m c — p x l d and m c w lGeV respectively. Both are in- 
teresting scales for physics in our world. Note the scale 
p x l d is associated with both the dark energy and the neu- 
trino massfrij]. One way m c could enter is through that 
fact (illustrated in Fig. [5]) that local Lorentz symme- 
try breaking theories are preferred for m very close to 
m c . We find it intriguing that other authors have associ- 
ated neutrino mass with Lorentz- violation [ll|, while our 
scheme appears to associate Lorentz-violation with the 
neutrino mass scale. 

The idea that physics at its most fundamental is sim- 
ply described by a random Hamiltonian is a curious idea 
which we motivated in previous work on the clock ambi- 
guity. There is much still to be explored before such an 
idea could look seriously viable. Still, we find it interest- 
ing that this idea has passed as many tests as it has so far. 
In particular, this paper shows that the random Hamilto- 
nian picture predicts that we would interpret the physi- 
cal world using Lorentz symmetric field theories. We hnd 
this result quite remarkable. 
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There, assumption 1 was dropped to include a sub- 
extensive Casimir energy component Ec while keep- 
ing the no-scale property (g — in assumption 2). 
This leads to the generalized Cardy formula S ~ 
V 1/d y/E c (E - E c ) . We could general ize further our ap- 
proach to include such a non-extensive component to ar- 
rive at the analog expression when the no-scale property 
is lost. 
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our energy spectrum test when an alternate form of the 
density of states was employed. This form was based on 
[l2T ] (as we mentioned briefly in section |4T2"|) . One would 
need to extend results of to include generalized dis- 
persion relations to fully explore the kv = (Vu) 1 ^ case. 



